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£Sj ■ We present a theoretical estimate for the cross-section of exclusive p + p and p°p°-meson production 

, ' in two photon collisions when one of the initial photons is highly virtual. We focus on the discussion 

of the twist 4 contributions which are related to the production of an exotic isospin 2 resonance of two 
p mesons. Our analysis shows that the recent experimental data obtained by the L3 Collaboration 
' at LEP can be understood as a signal for the existence of an exotic isotensor resonance with a mass 

\ around 1.5 GeV. 

>. 1 I. INTRODUCTION 

t> : 

, Exclusive reactions 7*7 — > A + B which may be accessed in e + e~ collisions have been shown 1\ to have a 
partonic interpretation in the kinematical region of large virtuality of one photon and of small center of mass 
energy. The scattering amplitude factorizes in a long distance dominated object - the generalized distribution 
amplitude (GDA) - and a short distance perturbatively calculable scattering matrix. A phenomenological 
analysis of the tttt channel Q has shown that precise experimental data could be collected at intense e + e~ 
collider experiments such as BABAR and BELLE. Meanwhile, first data on the p°p° channel at LEP have been 
Oh, published [3j and analyzed showing the compatibility of the QCD leading order analysis with experiment 
q; at quite modest values of Q 2 . 

In this paper, we focus on the comparison of processes 7*7 — * p°p° and 7*7 — * p + p~ in the context of 
1-^ ■ searching an exotic isospin 2 resonance decaying in two p mesons; such channels have recently been studied 
\ at LEP by the L3 collaboration 0,0- A related study for photoproduction [f| raised the problem of p°p° 
enhancement with respect to p + p~ at low energies. One of the solutions of this problem was based on the 
prediction Q and further exploration Q of the possible existence of isotensor state, whose interference with 
the isoscalar state is constructive for neutral mesons and destructive for charged ones. This option was also 
independently considered in [{| . The crucial property of such an exotic state is the absence of qq wave function 
at any momentum resolution. In other words, quark-antiquark component is absent both in its non-relativistic 
description and at the level of the light-cone distribution amplitude. This is by no means common: for instance, 
the 1 h state which is a quark-gluon hybrid at the non-relativistic level is described by a leading twist quark- 
antiquark distribution amplitude (lfjj . Contrary to that, an isotensor state on the light cone corresponds to the 
twist 4 or higher and its contribution is thus power suppressed at large Q 2 . This is supported by the mentioned 
L3 data, where the high Q 2 ratio two of the cross sections of charged and neutral mesons production points out 
an isoscalar state. 

We studied both perturbative and non-perturbative ingredients of QCD factorization for the description of an 
isotensor state. Namely, we calculated the twist 4 coefficient function and extracted the non-perturbative matrix 
elements from L3 data. Our analysis is compatible with the existence of an isotensor exotic meson with a mass 
around 1.5 GeV. 
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II. AMPLITUDE OF 7*7 -> pp PROCESS 

The reaction which we study here is e(k) + e(l) — > e(fc') + e(l') + p(j>±) + p(p2), where p stands for the triplet 
p mesons; the initial electron e(k) radiates a hard virtual photon with momentum q — k — k' , with q 2 = —Q 2 
quite large. This means that the scattered electron e(k') is tagged. To describe the given reaction, it is useful to 
consider the sub-process e(fc) + 7(5') — * e(k') + p{p\) + p{p2)- Regarding the other photon momentum q 1 = I — I', 
we assume that, firstly, its momentum is almost collinear to the electron momentum I and, secondly, that q 12 
is approximately equal to zero, which is a usual approximation when the second lepton is untagged. 
In two p meson production, we are interested in the channel where the resonance corresponds to the exotic 
isospin, i.e I — 2, and usual J PC quantum numbers. The J PC quantum numbers are not essential for our 
study. Because the isospin 2 has only a projection on the four quark correlators, the study of mesons with the 
isospin 2 can help to throw light upon the four quark states. We thus, together with the mentioned reactions, 
study the following processes: e(k) + e(l) — » e(fc') + e(V) + R(p) and e(k) + j(q') — > e(k') + R(p), where meson 
R(p) possesses isospin 7 = 2. 

Considering the amplitude of the 7*7 subprocess, we write 

%i)(77* - PP) = < W 4 J) / d 4 «i^e-*' / -* 1 - < «-* 9 <p(p 1 )p(p 2 )|r[J |t (« 1 )J I/ (« 2 )] |0), (1) 

where denotes the quark electromagnetic current J M = ipQl^ with the charge matrix Q belonging to 
SUf(2) group. The photon polarization vectors read 

■; w -(*&34 e —(^4 ^=(J?' 0A ^)' 

for the real and virtual photons, respectively. The coefficient functions of twist 2 operators to Operator Product 
Expansion of currents product in Q were discussed in detail in while the contributions of new twist 4 
operators are described by coefficient functions calculated long ago in when considering the problem of 
twist 4 corrections to Deep Inelastic Scattering. 

Let us now turn on the flavour or isospin structure of the corresponding amplitudes. The pp state with 1 = 
can be projected on both the two and four quark operators, while the state with I = 2 on the four quark 
operator only. Indeed, let us start from the consideration of the vacuum-to-pp matrix element in Q 

(pV 1^(0)1^)10) = S ab I fg ^-° + is abc r c fg ^=\ (3) 

where the quark fields are shown with free flavour indices and T stands for the corresponding 7-matrix. The 
isoscalar and isovector GDA's <& 7 in are well-known, see for instance [l2^. Note that, in ©, the correspon- 
dence between triplets {p 1 , p 2 , p 3 } and {p + , p~ , p } is given by the standard way. 

Moreover, for the coefficient function at higher order in the strong coupling constant, the corresponding matrix 
element gives us 

(p a p b I [4>h (o)r^ ffl (r))} (z)T 2 ^ g2 (0]|o> . (4) 



|/=o,/,=o (5) 



$'=2,/ 2 =0 (6) 

for the isospin 2 and I z = projection of the four quark operator in Q). The four quark GDA's $ / ' / ^ =0 can be 
defined in an analogous manner as the two quark GDA's. Hence, one can see that the amplitudes JQ) for p°p° 
and p + p~ productions can be written in the form of the decomposition: 

a _ ,7=0 : 7,=0 , .7=0,7^=0 , ,7=2,7,=0 /-x 

+ ) ~~ ■ H (+,+)2 ~ l ~" H (+,+)4 + " a (+ : + )4 ' \' ) 



Using the Clebsch-Gordan decomposition, we obtain 

_ 7=0 : 7 Z =0 
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for the isospin and I z = projection of the four quark operator in 0), and 
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where the subscripts 2 and 4 in the amplitudes imply that the given amplitudes are associated with the two and 
four quark correlators, respectively. The amplitudes corresponding to p + p~ production are not independent 
and can be expressed through the corresponding amplitudes of p°p° production. Indeed, one can derive the 
following relations: 

4:°+) 7 f °(77* - P + P-) = °(77* - P°P°) for k = 2, 4 

°(77* - p + p -) = -^;: 2 + ) / r°(77* - P o P o ). m 

The amplitude of two p meson production in two photon collision can be also presented through a resonant 
intermediate state. The vacuum to pp matrix element in the r.h.s. of can be traded for 

{p{Pi) P (P2)\R I (p)) 2 _ p2 V r (fl J (p)|T [^(0)^(^10). (9) 

7=0,1,2 R 1 P R ' RI 

where R 1 (p) is the resonance with three possible isospin 7 = 0, 1, 2. Note that, in our case, only isospin and 2 
cases are relevant due to the positive C-parity of the initial and final states. The matrix element (p p\R T ) defines 
the corresponding coupling constant of meson and (R T \T [J M (0) J^(z)] |0) is considered up to the second order 
of strong coupling constant as, i.e this matrix element is written as a sum of two- and four-quark correlators. 



III. DIFFERENTIAL CROSS SECTIONS 



Previously, the theoretical description of the experimental data collected for the p°p° production has been 
performed in Now, the subject of our study is the differential cross section corresponding to both the p°p° 
and p + p~ productions in the electron-positron collision. 

Using the equivalent photon approximation |l3| we find the expression for the corresponding cross section : 



d(7 ee — >eepp 

dQ 2 dW 2 



J ■■ J dcosO d(j)dx2—Fww{x2) 



da, 



dQ 2 dW 2 dcos6 < 



(10) 



where the usual Weizsacker- Williams function Fww is used. In (|10fl . the cross section for the subprocess reads 



da, 
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(11) 



where the amplitude Ar +<+ \ is defined by J7J. For the case of p°p° production, the cross section H1U|) takes the 
form: 



da ee ^ eep o p o _ 100a 
dQ 2 dW 2 ~ 9 

/ r R o Mfto 



G{S eei Q\W 2 )P 
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+ Q 2 S4 



as(Q 2 )M 2 R2 cI=2j3 ^ 

Q 2 4 



(12) 



where T R i stand for the total widths. The dimensionful structure constants S 



i,h=o 



and S 



7=0,7 3 =0 



are related 



to the nonperturbative vacuum-to-meson matrix elements. The /3-functions are also defines in the standard 
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ways: = J 1 - 4m 2 /W 2 and /3j = ^/l - Am 2 /M RI . The function G in (O is equal to 

l 

G(5 ee ,g 2 ,^ 2 ) = /" dx 2 F WW (x 2 ) 



x 2 S 2 e Q 2 x 2 S ee Q 2 (Q 2 + W 2 ) Q 2 {Q 2 + W 2 ) 2 

The differential cross section corresponding to p + p~ production can be obtained using (jSJ, we have 
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(13) 



(14) 



Note that we have explicitly separated out, in (|f 2|l and (|f 4|l . the running coupling constant as(Q 2 ) which 
appears in the twist 4 terms. Because of we will study the Q 2 dependence of the corresponding cross sections at 
rather small values of Q 2 , we use the Shirkov and Solovtsov's analytical approach 01 to determine the running 
coupling constant in the region of small Q 2 . Detailed discussion on different aspects of using the analytical 
running coupling constant may be found in [lflL flfij ] and references therein. 



IV. LEP DATA FITTING 



In the previous section we derived the differential cross sections da ee ^ eepp /dQ dW for both the p p and 
p + p~ channels, based on the QCD analysis. These expressions contain a number of unknown phenomenological 
parameters, which are intrinsically related to non perturbative quantities encoded in the generalized distribution 
amplitudes. One should now make a fit of these phenomenological parameters in order to get a good description 
of experimental data. The best values of the parameters can be found by the method of least squares, \ 2 ~ 
method, which flows from the maximum likelihood theorem, but we postpone a comprehensive % 2 -analysis to a 
forthcoming more detailed paper. Here, we implement a naive fitting analysis to get an acceptable agreement 
with the experimental data. Thus, we have the following set of parameters for fitting: 

p = {m r0 , r>, m r2 , r fl2 , s^= < /3 =°, sf°< /3=0 , sf 2j3=0 }. (is) 

We start with the study of the W dependence of the cross sections. For this goal, following the papers 
we determine the cross section of process ee — > eepp normalized by the integrated luminosity function: 



lT {(w)) = J - -* y r Z^,r/A, ■ (is) 



/ dQ 2 C{Q 2 ,{W))a^{Q 2 , (W)) 
JdQ 2 C(Q 2 ,(W)) 

where the definition of the luminosity function C is taken from [3. The value (W) corresponds to the center 
of each bin, see |3|, |5j ■ Focussing first on the region of larger Q 2 we fit the parameters associated with the 
dominant contribution which comes from the twist 2 term amplitude, which is associated with the non-exotic 
resonance (or background) with isospin 1 = 0. Generally speaking, there are many isoscalar resonances with 
masses in the region of 1 — 3 GeV. To include their total effect we introduce a mass and width for an " effective" 
isoscalar resonance. We then determine the values of the mass and width by fitting the data for the region of 
larger Q 2 (i.e., when Q 2 is in the interval 1.2 < Q 2 < 8.5 GeV 2 ). We thus can fix the parameters S2 =0 ' /3=0 , 
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Mro and Tflo. Good agreement can be achieved with M R o = 1.8 GeV, T R o = 1.00 GeV and S^ -0 ' 73-0 within 
the interval (0.12, 0.16). As can be expected the width of the effective isoscalar "resonance" is fairly large. It 
means that we actually deal with a non-resonant background. 

Next, we fit the W— dependence of the cross section for small values of Q 2 , i.e. 0.2 < Q 2 < 0.85 GeV 2 . In 
this region all twist contributions may be important. We find that the experimental data can be described by 
the following choice of the parameters: M R 2 = 1.5 GeV, T R 2 = 0.4 GeV while the parameters S^ -0 ' 73-0 and 
S4 =2 ' /3=0 are in the intervals (0.002, 0.006) and (0.012, 0.018), respectively. 




Figure 1: W dependence of the cross section <J ee ^ eep o p o normalized by the integrated luminosity function, in the 1.2 < 
Q 2 < 8.5 region. The short-dashed line corresponds to the leading twist 2 contribution; the dash-dotted line to the twist 
4 contribution; the middle-dashed line to the interference of twist 2 and 4 contributions. The solid line corresponds to 
the sum of all contributions. 




Figure 2: W dependence of the cross section cr ee _, eep + p - with the same conventions as for Fig. 

Further, we include in our analysis the Q 2 dependence of p°p° and p + p~ production cross sections, i.e. 
da ee ^ eep p j dQ 2 , which should fix the remaining arbitrariness of the parameters. We finally find that the best 
description of both W and Q 2 dependence is reached at 

M R 2 = 1.5 GeV, T R 2 = 0.4 GeV, 

S'f 0,73=0 = 0.12 GeV, S 7=0 ' 73=0 = 0.006 GeV, S 7=2 ' /3=0 = 0.018 GeV. (17) 

Note that these rather small values of twist 4 structure constants S4 compared to the twist 2 structure constant 
S2 indicate that leading twist contribution dominate for the values Q 2 around or greater than 1 GeV 2 . This 
should be compared with what was obtained in a particular renormalon model in [l7J • 

Our theoretical description of the LEP experimental data are presented on Figs. The plots depicted on 

Figs. EUHhave the following notations: the short-dashed line corresponds to the contribution coming from the 
leading twist term of l|12(l ; the dash-dotted line - to the contribution from the twist 4 term of ljl2|l ; the middle- 
dashed line - to the contributions from the interference between twist 2 and twist 4 terms of (|12H and l|14f> : the 



G 



long-dashed line - to the contribution from the interference between isoscalar and isotensor terms. Finally the 
solid line corresponds to the sum of all contributions. On Fig. [5] we present the LEP data and our theoretical 
curves for both the p°p° and p + p~ production differential cross sections as functions of Q 2 . The solid line on 
Fig. corresponds to the p°p° differential cross section while the dashed one - to the p + p~ differential cross 
section. 
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Figure 3: W dependence of the cross section a, 



, _ eep o p o normalized by the integrated luminosity function in the 0.2 < 
Q 2 < 0.85 region. The short-dashed line corresponds to the leading twist 2 contribution; the dash-dotted line to the 
twist 4 contribution; the middle-dashed line to the contributions from the interference between twist 2 and twist 4 terms; 
the long-dashed line to the contribution from the interference between isoscalar and isotensor terms. The solid line 
corresponds to the sum of all contributions. 
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Figure 4: same as Fig. I^Jfor the cross section cr ee ^ ee „+ p - in the 0.2 < Q 2 < 0.85 region. 



V. DISCUSSIONS AND CONCLUSIONS 



The fitting of LEP data based on the QCD factorization of the amplitude into a hard subprocess and a generalized 
distribution amplitude thus allows us to claim evidence of the existence of an isospin 1 = 2 exotic meson 0,1110] 
with a mass in the vicinity of 1.5 GeV and a width around 0.4 GeV. The contributions of such an exotic meson 
in the two p meson production cross sections (see, (fT2*)l and are directly associated with some twist 4 

terms that we have identified. At large Q 2 , these twist 4 contributions become negligible and the behaviours 
of the p°p° and p + p~ cross sections are controlled by the leading twist 2 contributions, see Fig. ^ and [21 Figs. 
[3] and 0] show the increasing role of higher twist contributions when decreasing Q 2 . Namely, the interference 
between twist 2 and 4 amplitudes gives the dominant contributions to p°p° production in the lower Q 2 interval, 
and is thus responsible of the W dependence of the cross section in these kinematics. In particular, in this 
interference term the main contribution arises from the interference between isoscalar and isotensor structures, 
see the long-dashed lines on Fig. |3]and0] 
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Figure 5: The Q 2 dependence of the differential cross sections da ee _ e> , p o p o / 'dQ 2 and dcr ee ^ eep + p - /dQ 2 . The solid line 
corresponds to the case of p°p° production; the dashed line to the case of p + p~ production. 

Analysing the Q 2 dependence, we can see that due to the presence of a twist 4 amplitude and its interference 
with the leading twist 2 component, the p°p° cross section at small Q 2 is a few times higher than the p + p~ 
cross section, see Fig. \5\ While for the region of large Q 2 where any higher twist effects are negligible the p°p° 
cross section is less than the p + p~ cross section by the factor 2, which is typical from an isosinglet channel (see 
also <P|) and ijTIjl ). 

The reaction 7*7 — > pp and its QCD analysis in the framework of Ref. thus proves its efficiency to reveal 
facts on hadronic physics which would remain quite difficult to explain in a quantitative way otherwise. The 
leading twist dominance is seen to persist down to values of Q 2 around 1 GeV 2 . Other aspects of QCD may 
be revealed in different kinematical regimes through the same reaction Its detailed experimental analysis 
at intense electron colliders within the BABAR and BELLE experiments is thus extremely promising. Data at 
higher energies in a future linear collider should also be foreseen. 

Note that the non-perturbative calculations of the relevant 1 = 2 twist 4 matrix elements also deserve special 
interest. In particular, one may follow the ideas developed for pion distribution |20| which allowed to relate 
higher and lower twists in multicolour QCD. The generalization for the case of p mesons, anticipated by the 
authors of l2fi . and use of crossing relations between various kinematical domains provided by Radon transform 
technique [21| may allow to apply these result in the case under consideration. 

In conclusion, let us stress that the L3 data allows to estimate the contribution of higher twist four quark light 
cone distribution to the production amplitude of vector meson pairs. Our numerical analysis leads to a rather 
small width for the corresponding resonant state, which is nothing else as an exotic four-quark isotensor meson. 
At the same time, a more elaborate experimental, theoretical and numerical analysis is required to confirm, 
with better accuracy, the smallness of the width and the existence of an exotic meson. 
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